Phantom-like behaviour in a brane-world model with curvature effects 
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Recent observational evidence seems to allow the possibility that our universe may currently be 
under a dark energy effect of a phantom nature. A suitable effective phantom fluid behaviour can 
emerge in brane cosmology; In particular, within the normal non self-accelerating DGP branch, 
without any exotic matter and due to curvature effects from induced gravity. The phantom-like 
behaviour is based in defining an effective energy density that grows as the brane expands. This 
effective description breaks down at some point in the past when the effective energy density becomes 
negative and the effective equation of state parameter blows up. In this paper we investigate if the 
phantom-like regime can be enlarged by the inclusion of a Gauss-Bonnet (GB) term into the bulk. 
The motivation is that such a GB component would model additional curvature effects on the brane 
setting. More precisely, our aim is to determine if the GB term, dominating and modifying the 
early behaviour of the brane universe, may eventually extend the regime of validity of the phantom 
mimicry on the brane. However, we show that the opposite occurs: the GB effect seems instead to 
induce a breakdown of the phantom-like behaviour at an even smaller redshift. 
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I. INTRODUCTION 

In the context of cosmology, one of the most relevant 
astronomical observations of the last decade are those 
from distant type la supernova implyingthat the uni- 
verse is in a state of accelerated expansion [l|, Q which has 
been latter on confirmed by other observational probes 
S 0, IE B 0, H| ■ The fundamental nature of what is driv- 
ing the cosmic acceleration is unknown, although many 
theoretical propositions have been put forward 0, llfjj ]. 
Invoking a cosmological constant to explain the late-time 
acceleration of the universe turns out to be the most eco- 
nomical option which moreover is in agreement with all 
the observational data Q. However, it turns out that the 
expected theoretical value of the cosmological constant is 
about 120 orders of magnitude larger than the measured 
one [3. 

Whatever it is the fuel inducing the late-time acceler- 
ation of our universe, from a phenomenological point of 
view and in the framework of general relativity, it can be 
described through a dark energy component (pd,Pd) with 
an effective equation of state w e s = Pd/ Pd- The current 
value of w e ff is extremely close to —1 Q; i.e. a cosmolog- 
ical constant equation of state, but it can be larger than 
-1 like in quintessence models or even smaller Q. 

The latter case disclosed above is of no lesser impor- 
tance. Quite on the contrary: If dark energy has a phan- 
tom nature, i.e., w e g < —1 [Tl| . then the scientific com- 
munity is facing a most considerable twofold challenge: 
(i) explaining the cause of the recent speed up of our uni- 
verse and (ii) also how to accommodate a phantom en- 



ergy component in our theoretical framework; i.e. what 
can cause w e s < — 1 without invoking a real phantom en- 
ergy which is known to violate the null energy condition 
and induce quantum instabilities 1 [13]. 

A possibility to mimic such a phantom-like behaviour 
is within the paradigm of string inspired brane-world 
models, where matter (standard model particles) is con- 
fined on a 4-dimensional (4D) hypcrsurfacc embedded in 
a higher dimensional space-time (the bulk) and where 
gra vity is the only interaction experiencing the full bulk 
[131 fl4| . More precisely, Sahni and Shtanov proposed 
the ADGP model with a phantom-like behaviour on the 
brane fl5l . Ilfl ] , without the need of any matter that vio- 
lates the null energy condition on the brane. 

The ADGP scenario is a 5D brane-world model with 
infrared modifications to general relati vity caused by an 
induced gravity term on the brane [U H [13, HO, 112 HH 
HI, [HI . The model is based on the normal or non self- 
accelerating branch of the Dvali-Gabadadze-Porrati pro- 
posal (DGP) [13, [3 j which, unlike the self-accelerating 
DGP branch is free from the ghost problem (25[, be- 
ing the brane filled with cold dark matter (CDM) and 
a cosmological constant which drives the late-time ac- 
celeration of the brane. The phantom-like behaviour is 
a consequence of the extra-dimension which screens the 
brane cosmological constant and it is based in mapping 
the brane evolution to that of an equivalent 4D general 
relativistic phantom energy model [l5l . fl~6j ] . More pre- 
cisely, the basis of this mimicry is an effective energy den- 
sity (in the 4D general relativistic picture) corresponding 
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1 We are referring here to a phantom energy component described 
through a minimally coupled scalar field with the wrong kinetic 
term. 
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to the cosmological constant corrected by the curvature 
effect due to the induced gravity term on the brane. This 
effective energy density grows as the brane expands and 
therefore effectively it behaves as a phantom fluid; i.e. 
w c ff < — 1, where w g corresponds to the ratio between 
the effective energy density and the effective pressure. 
The ADGP model in [HI, [l6[ is by far the simplest way to 
mimic a phantom- like behaviour in a brane- world setup 2 . 
Other brane proposals aiming to produce such a mimicry 
are based on a bulk filled with matter and/or on an en- 
ergy exchange between the brane and the bulk, therefore 
modifying the effective equation of state of dark energy 
on the brane [27| ■ 

The effective description of the phantom behaviour 
in the ADGP model breaks down at a finite redshift 3 
(cf. figure 6 and section IV for a more detailed descrip- 
tion); i.e. the effective energy density vanishes and be- 
comes negative over a certain redshift. When the effective 
energy density vanishes, the effective equation of state 
blows up. Given that the phantom-like behaviour results 
from (i) induced gravity effects on the brane causing cur- 
vature corrections and (ii) describing the brane model as 
a 4D relativistic phantom energy setup, could the break 
down of the phantom-like behaviour be eliminated by 
considering further curvature effects on the brane-world 
scenario? This is the main question we address in this 
paper. 

We will model such additional and new curvature ef- 
fects through a Gauss-Bonnet term (GB) in the bulk 
[30l . l3ll ] . The reason behind including this specific curva- 
ture terms is that it induces an ultraviolet correction on 
the brane [3(| HH , as expected from high-energy stringy 
features, and therefore may modify the phantom-like be- 
haviour at earlier times. Eventually, affecting its long 
term dynamics and even possibly preventing the men- 
tioned break down to occur at all. The other reason 
for considering such a curvature term was anticipated 
recently 0, [23] : even though the DGP model is char- 
acterised by an interesting infrared effect of gravity oc- 
curring with respect to general relativity, which for the 
self- accelerating branch can lead to a late-time accelera- 
tion on the brane even in the absence of any exotic matter 
invoked to produce the dark energy effect [l||], it would 
be expected that a consistent DGP brane model would 
have also ultraviolet modifications as well, associated to 
high-energy stringy effects at earlier times. 

This paper is therefore outlined as follows. In section 
HI! we define a brane-world model, henceforth designated 
as ADGP-GB model, and constrain the set of parame- 
ters that defined it in such away that the brane is cur- 
rently accelerating. We also comment on the non super- 



2 It is also possible to mimic a phantom-like behaviour in a con- 
sistent way in scalar-tensor dark energy models 26]. 

3 We would like to stress that the breakdown of the effective phan- 
tom description in the ADGP model does not imply any sort of 
singularity on the brane nor in the bulk. 



acceleration of the brane. In section IIIII we solve the 
cubic Fricdmann equation for the normal DGP branch 
with a GB term in the bulk: The reason is to obtain 
an accurate description of the effective energy density, 
that will behave like a phantom component on the brane, 
which depends explicitly on the Hubble rate. In section 
IIVI we subsequently show how a mimicry of a phantom 
behaviour takes place on the brane without considering 
any matter that violates the null energy condition on the 
brane. Then, we compare the behaviour of the ADGP- 
GB setting with the behaviour found on the ADGP model 
flU [H[ . Finally, in section [V] we summarise and con- 
clude. We also present some results related to the so- 
lutions of the cubic Friedmann equation (fT7|) in the ap- 
pendix A. On the other hand, in the appendix B, we 
show under which conditions the ADGP model is recov- 
ered from the model we propose. 



II. ACCELERATING ADGP-GB MODEL AND 
PARAMETER CONSTRAINTS 

The generalised Friedmann equation of a brane with 
induced gravity embedded in a 5D Minkowski bulk with 
a GB term reads 4 [H, H<| 

(l+^-aH^j H"=(r c H^-^p^\ (1) 

where a mirror symmetry has been assumed across the 
brane. In the previous equation r c is the crossover scale in 
the DGP model [Tt} and has length unit. This parameter 
measures the strength of the induced gravity effect on 
the brane and is related to the 4D and 5D gravitational 
constants by 



On the other hand, the parameter a measures the 
strength of the GB curvature effect on the brane and 
has length square unit and is positive [30I ]. 

If a = 0, then the induced gravity in the DGP setup 
modifies the late-time evolution of the brane with respect 
to the standard 4D relativistic case [l8| (for an alter- 
native approach where the induced gravity effect corre- 
sponds to a correction to RS model at high energies see, 
e.g., H^). However, if it is instead r c = then the GB 
term modifies the early-time evolution of the brane [30| . 

Equation JT]) can be conveniently rewritten as 

ff 2 = |p±l(l + ^ 2 )ij, (3) 



4 We restrict to the Friedmann equation that has an induced grav- 
ity limit and therefore contains the DGP model [29| . 
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which generalises the Friedmann equation of the self- 
accelerating DGP solution [13, Ell (+ sign in Eq. ^ 
with a = 0); i.c we recover the self-accelerating solution 
when o = 0. On the other hand, Eq. also has as a 
particular solution the DGP normal branch or non-self- 
accelerating solution [It], EH (— sign in Eq. (J3J) when 
a = 0). 

From now on, we restrict to the normal branch; i.e. — 
sign in Eq. (J3|). In addition, we consider that the energy 
density of the branc p corresponds to a CDM component 
with energy density p m and a cosmological constant A 



P = Pn 



A, 



(4) 



with the latter driving the late-time acceleration of the 
brane. We will refer to this scenario as the ADGP-GB 
model. The total energy density of the brane is conserved 
and therefore also the sector corresponding to the CDM, 
which scales in the standard way with the redshift 



p mo {l + zf 



(5) 



From now on a subscript stands for the observed cur- 
rent value of a given quantity. Finally, the Friedmann 
equation on the brane can be presented as 

e 2 {z) = o ro (i + zf + n A - 2^/aT c [i + n a E 2 ( z )] e(z), 

(6) 

where E(z) = H/Hq and 



2 

3H 2 



Oa = 



«|A 



(I) 



are the usual convenient dimensionless parameters while 
the new parameter £l a is defined as 



(8) 



Evaluating the Friedmann equation © at z = gives 
a constraint on the cosmological parameters of the model 



n A 



(9) 



For fl a — we recover the constraint in the ADGP 
model. The constraint (J9j> implies that the region f2 TO + 
f^A < 1 is unphysical. Moreover, although the brane 
is spatially flat, the previous constraint can be inter- 
preted in the sense that our model constitutes a mimic 
of a closed FLRW universe in the (£l m ,£l A ) plane. In 
particular, this is likewise to what happens in ADGP, 
QDGP and CDGP models 0, QJ, El, [13, H . We re- 
call that the QDGP and CDGP models correspond to 
variants of the ADGP scenario, where the late-time evo- 
lution of the universe is driven by a quiessence [2l| and a 
Chaplygin gas (24j, respectively, instead of a cosmologi- 
cal constant. Their dark energy effect is more dynamical 
and the phantom divide (or the w = — 1 line) crossing is 
possible in the QDGP and CDGP unlike in the ADGP 
model [2ll , [24| . We remind that the interest on modelling 



a mimicry of a phantom divide crossing is based on the 
possibility (backed by recent observational data) that the 
equation of state may have crossed the cosmological con- 
stant barrier (w = —1). 

Coming back to our model, if the dimensionless 
crossover energy density Q ra is the same in a ADGP 
model and in our model (which is not necessarily the 
case), then the similarities with a spatially closed uni- 
verse are made more significant from the GB effect, since 
il a > 0. Notice that this statement also applies to the 
variants of the ADGP brane mentioned previously and 
their generalisations by the GB effect, if the acceleration 
of both branes is driven by the same sort of dark energy. 
In fact, in this case the appropriately modified constraint 
(|9|) would read 



n 



DE 



(10) 



where JIde correspond to the current dimensionless en- 
ergy density of dark energy on the brane which can be for 
example modelled by a quiessence or a Chaplygin gas. 

Furthermore, by imposing that the universe is cur- 
rently accelerating; i.e. the deceleration parameter q = 
— (H/H 2 + 1) is currently negative, where 



'in 



^(1 + Q a ) 



(11) 



we obtain another constraint on the set of cosmological 
parameters Q m , Q rc and Q ai which reads 



3a„ < 2 



112) 



An example of the cosmological evolution of the deceler- 
ation parameter is given in Fig. [1] where it can be seen 
that the brane accelerates at late-time. 

On the other hand, the modified Raychaudhuri equa- 
tion follows easily from the Friedmann equation of the 
brane and the conservation of the brane energy density. 
It can be written as 



H 



n m (i + z) 3 E{ z ) 



2 E{z) + JVlTX 1 + Z& a E 2 {z)) ' 



(13) 



where a dot stands for the derivative respect to the cos- 
mic time. The key point of the previous equation is that 
the brane never super-accelerates; i.e. the Hubble rate 
decreases as the brane expands. Nevertheless, as we will 
show in section ITVl a phantom- like behaviour takes place 
at the brane: This occurs without including any matter 
that violates the null energy condition. The phantom-like 
behaviour is based in defining an effective energy density 
which corresponds to a balance between the cosmological 
constant and geometrical effects encoded on the Hubble 
rate evolution. Therefore, in order to get the evolution of 
the effective energy density with the redshift it is neces- 
sary to solve the cubic Friedmann equation of the normal 
branch (Eq. ((3]) with (-) sign). The solutions and an anal- 
ysis of the mentioned Friedmann equation is presented on 
the next section. 
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where 



FIG. 1: Plot of the deceleration parameter, 
q = -(H/H 2 + 1), versus the redshift. The set 

(n m ,fi A ,fir c ,n a ) = (0.26,0.7602, 1CT 4 , 0.01). As can 
be seen the brane accelerates at late-time when q gets 
negative. 
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"3 \ b -9 



1 - V(l - 36) ; 



(21) 
(22) 



for the analysis of the number of physical solutions of 
the modified Friedmann equation (TTT|) . If A/" is positive 
then there is a unique real solution. On the other hand, 
if Af is negative there are 3 real solutions. Finally, if Af 
vanishes, all roots are real and at least two are equal. 



Case 1: < b < 



III. NORMAL DGP BRANCH WITH 
GAUSS-BONNET EFFECT 

Let us herein solve analytically the Friedmann equa- 
tion J3| for the (-) sign for the reasons mentioned previ- 
ously 5 . Concerning this aim, it is convenient to introduce 
the dimensionlcss variables 



8 a , 



(14) 



P 
b 



32 K 2 n 2 

t^^itp = [n A + Mi + ^) 3 ] . (is) 

At r c 
8 a 
3^ 



(16) 



Then the Friedmann equation can be rewritten as 

H 3 + H 2 + bH -p = 0. (17) 

The number of real roots is determined by the sign of the 
discriminant function Af defined as 13511 



Af = Q 3 + R 2 , 



where Q and R read 



It is helpful to rewrite Af as 



2 P 



(18) 



M = t(p- Pi)(p~ Pz), 



(20) 



This is by far the most interesting physical case as we 
expect b to be small because it is proportional to f2 rc [see 
Eq. (fT6|) ] and the equivalent quantity in the ADGP sce- 
nario is relatively small [13, ; We do not expect that 
f2 rc in our model to be very different from that in the 
ADGP model. Furthermore, from the mimicry of our 
model regarding a closed FLRW universe [see Eq. 
and the constraint on the curvature of the universe; for 
example from the recent WMAP 5 years data in combi- 
nation with the baryon acoustic oscillations Q , and 
f2 a should be small and therefore b is also expected to be 
small. 

The analysis of this case is slightly involved. The rea- 
son is essentially that if < b < j then pi and pi are real 
[see Eqs ([2"Tj) and (f2"2"j) ] and therefore it is not as straight- 
forward as in the next cases to know the number of real 
solutions of the cubic Friedmann equation on H ; i.e. to 
know the sign of M [see Eq. (|20|)] . 

More precisely, in this case p2 < and < p\ [see 
Eqs [|2"I ]) - (J2"2]) ]. Then, the number of real roots of the 
cubic Friedmann equation (|17[) depends crucially on the 
minimum energy density of the brane: 



(23) 



i.e. the asymptotic value of the total energy density at 
z = — 1. We enumerate next the possible different situa- 
tions [see Fig. [3]: 



1- pi < pmin 

The minimum energy density of the brane is such that 
Pi < Pmin- Then the function Af is positive and there is 
a unique solution. The condition pi < p m i n implies 



5 The Friedmann equation JlJ has been previously analysed in l28l . 
Our new analytical analysis is exclusively based on the conve- 
niently modified cubic Friedmann equation (|3j which make the 
study much easier because the discriminant of the cubic equa- 
tion, J\f defined in Eq. 1181 . is much simpler. 



1 

37^ 



2 r 



l + \/(l-3&)= 



(24) 



and therefore constrains the set of allowed values of Q rc , 
fl a and f^A- In Fig. [5] we show the set of those parameters 
that fulfil the inequality as the uncolourcd area while 
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b = 4fi r Sl a 



FIG. 2: The coloured area corresponds to the set (Or c , 
Oc^a) that does not fulfil the inequality (|24[) . This is the 
most likely situation as Q rc and Q, a are expected to be small. 
On the other hand, the uncoloured area correspond to the set 
(O rc ,n a ,f2A) that fulfil the inequality (|24|l . 



the red coloured area corresponds to the set (O ro , Sl a , 
SI a) that does not fulfil the condition ([24]) . 

Finally, the expansion of the brane is described by 
Eq. {SI) or 



Hi 



where 77 is defined as 



cosh(r/) 



i? 



sinh(r/) 



'Q 3 + R 2 



(25) 



(26) 



and 77 m i n < T). The parameter rj m i n is defined as in 
Eq. (|26p with p = p m j n and this value of 77 is reached 
at z = — 1. It turns out that the expanding brane so- 
lution is asymptotically de Sitter in the future. On the 
other hand, at early time (large 77) matter on the brane 
is dominated by dust though its cosmological evolution 

does not correspond to the standard relativistic dust case 

1 

because at high redshift H ~ p~m, where p m is defined as 
in Eq. (|15jl . This is a consequence of the dominance of 
GB effects at high energy. This feature applies also to 
the high energy regime described in the next subsection. 



2. pmin < pi 

The minimum energy density of the brane is such that 
Pmin < pi ■ Consequently, the inequality (|24|) is not satis- 
fied and this again restricts the set {Sl rc , O a , CI a) which 
in this case corresponds to the coloured area in Fig. [2j As 
this figure highlights this is the most likely situation as 
we expect Q ra and tt a to be small for the reasons stated 
before. 

As the energy density blue-shifts backward in times; 
i.e. the energy density grows backward in times we can 
distinguish three regimes: 



• Limiting regime: p — p\. 

• Low energy regime: p min < p < pi- 

During the high energy regime, the energy density of 
the brane, p, is bounded from below by p~\ and there- 
fore the function Af is positive or equivalently there is 
a unique solution of the cubic Friedmann equation (1 1 7|) . 
During this regime, the expansion of the brane is de- 
scribed by Eq. (|25|) where < 77 and defined in Eq. (|26j) . 
When 77 — ^ 0, the energy density of the brane approaches 

Pi- 

During the limiting regime, p = p\. Consequently M 
vanishes and there are two solutions: 



1 



H x = - ( 2VT^3b - 1 ) , 



H 2 = -- (VT^35+1 



(27) 
(28) 



The high energy regime connects with the limiting regime 
through Hi. The negative solution H2 is not relevant 
physically. 

Finally, at the low energy regime the total energy den- 
sity of the brane is bounded from above by p\ . Then M 
is negative and there are 3 different solutions [see Fig. [3]. 
One of this solution corresponds to an expanding brane 
while the other two corresponds to contracting branes: 

It can be shown that the expanding solution (H > 0) is 
described by Eq. (|A1[) and more appropriately rewritten 
as 



H, = 



where 



1 



2 vl— 35 cos 



- 1 



< 9 < 



cos(0) 



R 



sin(0) 




(29) 



(30) 



For 9^0, the energy density p approaches p\\ i.e. the 
low energy regime is connected with the high energy 
regime through the solution (|27|) . On the other hand, 
$max is defined as in Eq. (|30|) with p = p m i n where the 
brane reaches its asymptotic de Sitter regime at z = — 1 . 
Notice that as matter redshifts on the brane, the angle 
9 gets larger. On the other hand, in this model, if the 
cosmological constant vanishes then the maximum angle 
9 is given by 9q where 



, 9 M 1 

cos — = , 

3 J 2 V / T^3F 



(31) 



High energy regime: p\ < p. 



and therefore, the Hubble rate vanishes. This feature 
signals that this brane solution does not corresponds to 
a self-accelerating brane. 

For completeness, we write down the remaining two 
solutions of the Friedmann equation (fTT|) when < b < 
1/4 and p < p~\. As it was anticipated before, these 
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solutions describe contracting branes and correspond to 
the solutions given in Eqs. (|A2|) and (|A3jl . They read 



H, 



2\/l - 36 cos ( | 



2Vl - 36 cos f ^4-^ I 



, < 8 < 8 max , 
(32) 

, < 8 < 8 max , 
(33) 



respectively. Unlike the solution Hi , these two solutions 
are contracting because H2 and H3 are negative. It can 
be shown that H 2 < H 3 . Both solutions approaches the 
same Hubble rate in the past at 8 = corresponding to 
the limiting solution (f2"8]) . Finally, if there is no cosmo- 
logical constant on the brane, then in the far future (at 
z = —1) the angle 9 is given by Eq. (|31|). where H2 and 
H3 approach constant negative values. 

Before ending we would like to point out that it is 
only the expanding branch with Hubble rate Hi that has 
a phantom-like behaviour which we will describe in the 
next section. An example of the three different solutions 
of the Friedmann equation (fT7f for < b < j and p min < 
pi can be seen in Fig. © . 

Hi 




FIG. 3: Plot of the dimensionless Hubble rates Hi, H2 and 
S3 against the dimensionless energy density p. The red curve 
corresponds to Hi. The blue curve corresponds to Hi and 
the other one to Hz. As can be seen it is only the solution 
corresponding to Hi that can exist in the "far" past and it 
is also the only expanding solution. If the minimum energy 
density of the brane p m m is larger than pi only the solution 
Hi exits. In the opposite case; i.e. p m j n < pi, the three 
solutions exist. 



B. 



Case 2: \ < b < k 

4 — 3 



If 6 i < b < i then pi < and pi < which implies 
that Af is positive because the total energy density of 
the brane p is positive. Consequently, there is a unique 
real solution and the cosmological evolution of the brane 



is unique. The dimensionless Hubble parameter is given 
in Eqs. (f25j) -([26 |) and r) min is defined as in Eq. (|26|) with 
P = Pmin- When 77 approaches its minimum value, the 
Hubble rate is constant and positive; i.e. the brane is 
asymptotically de Sitter and its expansion is dominated 
by the cosmological constant. If all matter on the brane 
redshifts and the total energy density on the brane van- 
ishes at z = — 1 then at rj = ryo, where 



cosh 



\ 3 



1 



2\/T — 3T 



(34) 



the Hubble rate vanishes. This feature is in agreement 
with the fact that this solution does not correspond to a 
self-accelerating branch. 

At high energy (large 77) matter on the brane is domi- 
nated by dust. However, its cosmological evolution does 
not correspond to the standard relativistic dust case be- 



cause at high redshift H 
high energy. 



Pm due to the GB effects at 



C. Case 3: b = | 



This constitutes a marginal case where b = 1/3; i.e. 
Q = l/(120 r J. The modified Friedmann Eq. $TTJ) has 
a unique real solution because N > 0. This can be no- 
ticed easily by realising that pi = p 2 < when 6=1/3 
and therefore the right hand side of Eq. (|2"0|) is always 
positive. The dimensionless Hubble parameter given in 
Eq. (|A1|) can be rewritten in a simple way as 



H = l 
3 



(1 + 27/fP - 1 



(35) 



The brane is therefore asymptotically de Sitter in the fu- 
ture (p — > p m i n ). At high energy /earlier time, the matter 

on the brane is dominated by dust although the dimcn- 

1 

sionless Hubble parameter redshift as pm- 



D. 



Case 4: \ < b 



In this case, Eq. (jTTJ) has a unique real solution because 
Af > as pi and p~2 are complex conjugates when \ < b 
and therefore the right hand side of Eq. (|20j) is always 
positive. 

The dimensionless Hubble parameter is given in 
Eq. (|A1|) and can be rewritten as 



2V36^ 
where now rj fulfils 

cosh(?/) 



L ( I 



Vmin < ?/, (36) 



/ R 2 

1 + —^T, sinh(r?J 



R 



(37) 



The parameter pi vanishes at b = 1/4. 



and the parameter r) mm is defined as in Eq. (|37[) with 
P = Pmin- The brane is therefore asymptotically de Sitter 



7 



(z = —1 corresponds to r/ approaching r/ m ; n ) though there 
is no self-accelerating solution. 

For the sake of completeness, we point out that if all 
matter on the brane redshifts and the total energy density 
on the brane vanishes at z = —1 then at this rcdshift 
V = Vo, where the Hubble rate vanishes, in agreement 
with the fact that this solution does not correspond to a 
self-accelerating branch, and where 770 satisfies 



sinh ( 

* > 



1 



(38) 



Finally, we have that at high energy (large 77) the brane 

is dust dominated although its cosmological evolution 

does not correspond to the standard relativistic dust case 
_ 1 

because H ~ pf n where p m is defined as in Eq. (fT5|) . 

The asymptotic de Sitter regime of the brane in all 
the cases numerated is due to the presence of a cosmo- 
logical constant on the 4D hypersurface (unlike the self- 
accelerating solutions In addition, the fact that the 
brane is asymptotically dc Sitter implies that there is no 
big rip singularity in the future despite that the brane 
has a phantom-like behaviour at late-time as we show in 
the next section. 

Before proceeding into discussing how a phantom- 
like behaviour takes place on the brane framework we 
use, there is a point worthy to emphasise. The Gauss- 
Bonnet parameter can, in principle, have an arbitrary 
value. However, being also considered as a perturbative 
term arising from string theory, it is sensible that within 
the model discussed here, realistic cosmological solutions 
should coincide with the ADGP cosmology in the limit 
a — > 0. We therefore present in the appendix B how this 
is the case for one of the herein found solutions (Eq. (|29j)) 
and under which conditions the ADGP model is recov- 
ered from the model we propose. This procedure follows 
the analysis in Refs. [28|, [3l|. Moreover, it will be the 
solution (|29| that will be employed in the next section. 



IV. PHANTOM-LIKE BEHAVIOUR ON THE 
BRANE AND GAUSS-BONNET EFFECT 



The phantom-like behaviour on the brane is based in 
defining a corresponding effective energy density p e g and 
an effective equation of state with parameter w e g. More 
precisely, the effective description is inspired in writing 
down the modified Friedmann equation of the brane as 
the usual relativistic Friedmann equation so that 



H 2 



Trip™ + Pes); 



(39) 



i.e. to map the brane evolution in Eq. ^ to the equiv- 
alent 4D general relativistic phantom cosmology with 
Friedmann equation (|39[) . In the previous equation the 



0.6 
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FIG. 4: Plot of the dimensionless effective energy den- 
sity versus the redshift. The set (Q m , &a, f2r , = 
(0.26,0.7602, 10" 4 ,0.01) which corresponds to b = 4 x 10~ 6 
where b is defined in Eq. (|16[1 . The chosen values for 
(fi m , fiA, fir c ) are in agreement with the best observational 
fit of the ADGP model [231 ] and suitable for our model as we 
expect GB effects in our model to correspond to small correc- 
tions to the ADGP scenario. 



effective energy density p c g reads 



PeS 



= A- 



I + -aH 2 H, 



3H 2 



n A -2y^~(l + n a E 2 (z))E(z) .(40) 



This effective energy density corresponds to a balance 
between the cosmological constant and geometrical ef- 
fects encoded on the Hubble parameter. On the other 
hand, gravity leakage at late-time screens the cosmolog- 
ical constant like in the ADGP scenario [HI, [l6| ■ This 
phantom-like behaviour is obtained without any mat- 
ter violating the null energy condition and without any 
super-acceleration on the brane. We stress that the de- 
pendence of poff on the redshift is known analytically by 
means of the different solutions H of the cubic Friedmann 
equation presented on the previous section. 

In Fig. [4] we show an example of the evolution of the 
dimensionless effective energy density K 2 p e g/(3H 2 ). In 
this example b = 4 x 10 -6 . Our choice for the val- 
ues of b is based on the fact that observationally the 
most favourable set of solutions Hi are those such that 
< b < j and p m i n < p\ (we refer the reader to the previ- 
ous section). Furthermore, the chosen value b is obtained 
for Q r<: = 10 -4 which is in agreement with the best fit of 
the ADGP model (23| and suitable for our model as we 
expect GB effects in our model to correspond to small 
corrections to the ADGP setup. 

As in the ADGP model, it is possible to define an ef- 
fective equation of state or parameter w e g associated to 
the effective energy density as 



Pcff + 3-ff (1 + w c g )p cff = 0. 



(41) 



This effective equation of state is defined in analogy with 
the standard relativistic case. Then using Eq. (|4"0")) . we 
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FIG. 5: Plot of the effective equation of state versus the red- 
shift. The set (fi m , fi A , Or c , 0») = (0.26, 0.7602, 10" 4 , 0.01) 
the same considered in Fig. [4] 



obtain 
1 + w cS 



1 H(l + 8aH 2 ) 

Hp ef [ 



2 



(42) 



m a E 2 ( z )) 



3 e(z) [n A 



Q a E 2 (z))E(z)] 



Because the brane never super-accelerates, i.e. H < 
[see Eq. (fl~3"|) ]. we can then conclude that p e ff mimics the 
behaviour of a phantom energy component on the branc: 
I.e. 1 + w c s < 0, as long as the effective energy density 
p e g is positive [see Figs. |4] and [5]. 

We would like to point out that a phantom energy 
component can be defined in two ways: (i) any matter 
such that its equation of state fulfils p/p = w < — 1 
or (ii) any matter whose energy density grows when the 
universe expands. Both definitions arc equivalent as long 
as the universe expands and the energy density of the 
phantom component is positive. Here, we are assuming 
that the phantom energy is defined as in (i) and it turns 
out that the definition (ii) is automatically satisfied. In 
this model, the opposite situation does not hold, in other 
words if condition (ii) is satisfied (which is always the 
case in this model) it does not imply that condition (i) 
is fulfilled. The reason behind it is that p g is always 
positive, 

3 -H(l + 8aH 2 ), (43) 



Pefi 



K 2 r c 



and therefore p e g always grows as the brane expands in- 
dependently of its sign. However, 1 + w e fi changes its 
sign, although in an abrupt way, when p e g vanishes and 
becomes negative [see Eq. (T12")) and Fig. [5]. 

When the effective energy p e g vanishes, the mimicry 
of the phantom behaviour breaks down; i.e. the map- 
ping between the ADGP-GB model and the 4D relativis- 
tic phantom cosmology model with Fricdmann Eq. (|39|) 
is no longer valid , although the ADGP-GB model is 
well defined at any redshift. At the redshift z b , where 
Pes(zb) = 0, the Hubble rate is constrained to fulfil [see 
Eq. gO])] 

[i + n a E 2 (z b )] E(z b )-n A = o. (u) 



On the other hand, the Fricdmann equation (|39p implies 

E 2 (z b ) = Cl m (l + z b f. (45) 

By combining the last two equations, we find that the 
redshift at which the mimicry of a phantom behaviour 
breaks down reads [35l 



Zb 



— S- 



2 yj ^7' c 



(46) 



where 



S± 



(47) 

When the dimensionless energy density, f2 Q , associated 
to the GB geometrical effects is much smaller than fi rc , 
the redshift z b fulfils 

2^A c (l + z b )^ = n A -±nl^ + 2 (J^Lj . (48) 

Therefore, in the limiting situation Q a = 0, we recover 
that the phantom-like description breaks down also in 
the ADGP setup at same point in the past 7 : 



2^n~n~(i 



Zb) 



= n, 



(49) 



In fact, the effective phantom-like description breaks 
down in the ADGP scenario when the analogues effec- 
tive energy density vanishes and therefore the effective 
equation of state parameter blows up [33] [see Eq. (|42|) 
with Q a = 0] . This again points to the fact the mapping 
between the ADGP model and the 4D relativistic phan- 
tom cosmology model with Friedmann expansion (|39|) 
(for a = 0) is no longer valid, although the ADGP brane 
description remains valid 8 . 

The above behaviour raises the following possibility: 
Can the phantom behaviour break down in our model at 
a redshift z b (|46p . such that z b > z b l Our motivation 
is that such a GB component would model additional 
curvature effects on the brane setting. The subsequent 
aim is then to determine if the GB term, even if dom- 
inating and modifying the early behaviour of the brane 
universe, may eventually extend its effect towards later 
times. Thus, could the regime of validity of the phantom 
mimicry in the ADGP setting be extended by considering 
further curvature effects on the brane- world scenario? 

The answer will depend on the cosmological parame- 
ters that characterise both models, three in the ADGP 
scenario namely (£lmft A ,Q? c ) and four in the model we 



7 We use a tilde to define quantities in the ADGP model and hence 
distinguish them from the ones used in our model. 

8 We thank Y. Shtanov for pointing out this to us. 
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FIG. 6: In this plot we show how the phantom-like behaviour 
of a ADGP model is modified by the effect of a GB curvature 
term in the bulk. The effect is schematically shown by the red 
arrow. Hence, the GB effect translates into a sooner break- 
down of the effective phantom-like picture; i.e. Zb < Zb- The 
redshifts values Zb and Zb correspond to the moments when 
the effective energy density vanishes and becomes negative in 
the pure ADGP and the modified one with GB effect in the 
bulk, respectively. 



are analysing namely (Q m ^Qf i _,fl rc ,Q a ). In order to be 
able to give a definite answer, we will first assume three 
different cosmological situations. More general situation 
will be analysed by means of three dimensional plots (see 
Fig. [3. 



A. Fixed f2 m and Qa 

We first assume the same amount of dark matter fi m 
and dark energy Q A in the ADGP model and its twin 
with GB effect. Therefore, the difference between the 
ADGP scenario and ADGP-GB in this case is encoded on 
the geometrical effects quantified through the crossover 
scale and the GB parameter. By using the constraint 
equation (JS]), it turns out that the dimensionlcss energy 
density associated to the crossover scale cannot be the 
same in both models. In fact, the dimensionless energy 
density related to the crossover scale Q,? c m the ADGP 
model is related to fl rc in the ADGP-GB model by 



n? e = n rc (i + n a ) 2 



(50) 



Therefore, in this situation the crossover scale in the 
ADGP model will be larger than in the ADGP scenario 
with a GB term in the bulk. Let z b and Zb be the red- 
shift at which the effective energy density vanishes in the 
ADGP setup and in our scenario, respectively. Then, as 
the effective energy density vanishes at those redshifts 
Eq.(gn]) implies 

n A = 2^/nT c [i + n a E 2 (z b )]E(z b ), 

= 2^n~E(h), (51) 

which by using Eq. ([50]) translates into 



E(S b ) = l + n a E 2 (z b ) 

E(zb) i + n a 



(52) 



Now we recall that 1 < E(z b ) as the branc do not super- 
accelerate and < Zb- Therefore, the right hand side 
of the previous equation is larger than one, which im- 
plies that z b < z b because of Eq. (|4"5"|) and the analogous 
relation 



E(zb) = n A (i + z b ) 



(53) 



in the ADGP model. 

In conclusion, the phantom-like behaviour of the 
ADGP model with GB effect breaks down at a smaller 
redshift than the phantom-like behaviour on the ADGP 
model under the assumption that Q m and J7a are the 
same in the ADGP setup with and without GB effect in 
the bulk [see Fig. [6]. 



B. Fixed Qa and Q r<: 

We next assume the same amount of dark energy Ha 
and the same dimensionless energy density £l rr in both 
models. Then, the constraint equation |9]) implies that 
the amount of dark matter is slightly larger in the ADGP- 
GB brane than in the ADGP brane because of the GB 
effects 



(54) 



At the redshifts z b and z b the effective energy densities 
in the ADGP and the ADGP-GB models vanish, respec- 
tively, and Eq. (|40|) implies 



E{z h ) 
E(z b ) 



= i + n a E 2 { Zb ) > 1, 



(55) 



because we have assumed the same fl\ and O rc in both 
models. Now, Eqs @5J), ([521) an <i imply 



Zb 



1 



Zb 



> 



(56) 



Therefore, z b > z b because the amount of dark matter is 
larger in the ADGP-GB brane than in the ADGP brane 
[see Eq. (|54j)] . In conclusion, we have that the phantom- 
like behaviour breaks down sooner in the ADGP-GB 
model than in the ADGP setup [see Fig. [6]. 



C. Fixed Q m and Sl rc 

In this case we assume that both models have the same 
amount of dark matter D, m and the same dimensionless 
energy density f2 rc . Now, we use Eq. © to constrain the 
remaining free parameters of both models 



VL\ — f2« = 2a 



(57) 



The previous constraint implies that the amount of dark 
energy would be slightly larger in the ADGP-GB model 
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FIG. 7: Plot of the redshift Zt, at which the phantom-like behaviour in the ADGP-GB model breaks down. When f2 Q = 0, we 
obtain the redshift Zb at which the phantom-like behaviour in the ADGP scenario breaks down. As can be noticed Zb < Zb- 



than in the ADGP model in this case. Because the effec- 
tive energy density p c g vanishes at Zb and Zb, it implies 
a constraint on how different can be the amount of dark 
energy in both models [sec Eq. (|4"D1 ] 



E{z b ) - E(z b ) + 2^n r n a E 6 {z b ). 



(58) 



By combining Eqs. (fST)) . (|55| and recalling that z& > 0, 
we can conclude that 



E(z b ) - E(z b ) = fl a [E 3 {z b ) - 1] > 0; 



(59) 



i.e. the dimensionlcss Hubble rate E{zb) is larger than 
the dimensionlcss Hubble rate E(zb). Therefore, Zb < Zb 
[see Eqs ([35]) and l[5"5|)]; i.e. we reach a similar conclu- 
sion to the one presented in the previous subsections [see 
Fig. [6]. 

So, far we have compared the redshifts Zb, Zb where the 
mimicry of the phantom behaviour breaks down in the 
ADGP-GB and ADGP models, respectively, under the 
assumption that two of the parameters that characterise 
both models are equals. We relax this condition in the 
next subsection. 



D. More general situations 



In order to compare 
z b = z b (0 A , flm , flf c ) and 



further the redshifts 

Zb — Zb (f^A j j ^r c : ) 

at which the phantom-like behaviour breaks down in the 
ADGP and ADGP-GB models, respectively, we relax 
the conditions imposed in the previous subsection. 

We consider that only one of the parameters (^a, ^m, 
il rc ) is fixed in the ADGP-GB branc. We further assume 
that the given value of the fixed parameter is very close to 
that obtained by constraining the ADGP model with ob- 
servational data (H(z), CMB shift parameter and SNIa 
data) 9 H| 



0.26 ±0.05, a c <0.05. 



(60) 



Our second assumption is based on the fact that we are 
considering the GB effects as small corrections to the 
ADGP model. Then, once we have fixed the value of one 
of the parameters (Qa, Qm, ^r c ), there remain only three 
free parameters in the ADGP-GB model. One of these 
parameters can be fixed by means of the cosmological 
constraint equation ^ . Therefore, at the end we are left 
with only two free parameters and as we know analyti- 
cally Zb [see Eq. (|46|) ] we can do a three dimensional plot 
of Zb in terms of the two free parameters. Our results are 
shown in Fig. [7] 

On the left hand side plot in Fig. [7j we have chosen 
ft A = 0.74, its best fit value for the ADGP model [H[. 
Then for a given set (Q m ,S} a ), n rc is fixed by means 
of the constraint |9]). In this plot z b is retrieved when 
il a = 0. As can be seen for a given value fi ro the largest 
Zb corresponds to Q a = 0; i.e. Zb < Zb- 

On the middle plot in Fig. [7J we have followed a similar 
procedure. We have chosen O m = 0.26, its best fit value 
for the ADGP model [H[. Then for a given set (fi A , 
£l ra is obtained by means of the constraint f9j . In this 
plot Zb is retrieved when f2 Q = 0. As can be seen for a 
given value fl\ the largest Zb corresponds to Q a = 0; i.e. 
z b < z b - 

Finally, on the right hand side plot in Fig. [71 we have 
chosen tt rc = 10~ 4 which fulfils the constraint ((60)) . For a 
given set (f^A, fi a ), ^ m is fixed by means of the constraint 
(|9|). In this plot z b is retrieved when Q a = 0. It is clear 
form the plot that for a given value ^a the largest Zb 
corresponds to fl a = 0; i.e. z b < Zb- 

Hence, we have that in all the situations A-D anal- 
ysed above the phantom-like behaviour in the ADGP- 



As far as we know, the DGP-GB model has been constrained 



by cosmological observations in Ref. [33}. However, it turns out 
that the mentioned paper deals with the self-accelerating branch 
of the model which is different from the one we explore in our 
paper (the one that contains the ADGP-GB model). The model 
analysed in Ref. 33] reduces to the self-accelerating DGP branch 
when a — » while our model reduces to the ADGP model which 
is contained in the normal DGP branch. Therefore, as we are 
considering the GB term as a perturbation of the ADGP model, 
we take as a good approximation to consider the cosmological 
parameters of the ADGP-GB model close to those of the original 
ADGP model. 
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GB brane breaks down sooner (smaller redshift) than in 
the ADGP brane. In more precise terms, the phantom- 
like behaviour in the ADGP is regular in the interval 
[0, Zb) but breaks down for higher redshifts. On the 
other hand, the phantom-like behaviour in the ADGP- 
GB scenario works in the interval [0, Zb) but fails for 
higher redshifts. Because Zb < Zb, the GB effect actually 
makes smaller the interval [0, Zb) and therefore it does 
not help to improve the situation. Please note (i) we are 
analysing the phantom-like behaviour through the map- 
ping between the ADGP / ADGP-GB models and the 4D 
relativistic phantom cosmology models with Friedmann 
Eq. (|39| from today (z = 0) till Sb/zb and (ii) the break 
down of the phantom-like behaviour does not imply any 
sort of singularity in the brane-world models. 



V. SUMMARY AND CONCLUSION 

In this paper we have analysed in some detail the 
ADGP-GB model which corresponds to a 5D brane-world 
model where the bulk is a 5D Minkowski space-time. The 
model contains a GB term in the bulk I30l . l3l| and an 
induced gravity term on the brane (l9l . l28j |. Our analy- 
sis was performed for the normal or non self-accelerating 
branch which we have assumed to be filled by CDM and 
a cosmological constant, the latter driving the late-time 
acceleration of the brane. We have shown how the brane 
accelerates at late-times (cf. Eq. (Tj2j) and Fig.Q]). 

The attractive and promising feature of this model is 
the role of the extra-dimension: It induces a mimicry 
of a phantom behaviour without resorting to any mat- 
ter that violates the null energy condition on the brane. 
This phantom-like behaviour happens without any super- 
acceleration of the brane (see Eq. (fT3|) ) and, therefore, 
the brane does not hit a big rip singularity in its fu- 
ture. Indeed, the brane is asymptotically de Sitter. This 
regime is reached when the cosmological constant domi- 
nates over the CDM component. Our model reduces to 
the ADGP scenario [HI, [l6| if the GB corrections in the 
bulk are put aside. 

Our motivation for considering a GB correction to the 
ADGP is threefold: (i) it is known that the mimicry of a 
phantom behaviour on the ADGP model breaks down at 
some point in the past [32| . This happens when the map- 
ping between the ADGP and the 4D relativistic phantom 
cosmology model breaks down. More precisely, when the 
effective energy density that mimics the phantom-like be- 
haviour gets negative and therefore the corresponding ef- 
fective equation of state parameter w e g blows up (see 
Eq. (|42l) for a = where a is the GB parameter), (ii) 
GB effects induce ultraviolet corrections on the brane 
j30| . as expected from high-energy stringy features, and 
therefore may modify the phantom-like behaviour at ear- 
lier times and may alleviate the shortcome mentioned in 
the previous point. And [Hi), even though the DGP 
model is characterised by an interesting infrared effect of 
gravity occurring with respect to general relativity, which 



for the self-accelerating branch can lead to a late-time ac- 
celeration on the brane even in the absence of any exotic 
matter invoked to produce the dark energy effect [l8j], it 
would be expected that a consistent DGP brane model 
would have ultraviolet modifications as well, associated 
to high-energy stringy effects at earlier times [28j . 

The phantom-like behaviour on the brane is based on 
(i) writing down the modified Friedmann equation of the 
brane as the standard relativistic Friedmann equation 
(see Eq. ([55)) ) and (ii) defining a corresponding effec- 
tive energy density, p e s, which grows as the brane ex- 
pands, and an effective equation of state (cf. Eqs. (|4"01) 
and (|42|) and Figs. [4] and [5]). The effective energy den- 
sity corresponds to a balance between the cosmological 
constant and geometrical effects encoded on the Hubble 
rate. This was done by generalising the way a phantom- 
like behaviour is obtained in the ADGP model [HI, [l6| 
and other of its variants [U [24[ where dark energy has 
a more dynamical character. As p e s depends explicitly 
on the Hubble rate, in order to get its evolution with 
the redshift it is necessary to solve the cubic Friedmann 
equation (|6|). This was done in section ITTT1 This analysis 
has also allowed us to constraint the set of most likely 
values of the model (see Fig. [2]) and therefore to pick 
up the suitable cubic solution of the Friedmann equation 

It turns out that the phantom-like behaviour also 
breaks down in the ADGP-GB model, namely when 
Pcft = 0; i.e. when the brane cosmological constant 
balances the geometrical effects described in terms of 
the Hubble rate (see Eq. (|30])). This feature highlights 
that the mapping between the ADGP-GB model and 
the 4D relativistic phantom cosmology model ceases to 
be valid although the brane description remains valid. 
The redshift, Zb, at which this event happens is given in 
Eq. (|46|) and depends on the set (Q m ,QA,Q ra ,Q a ); i.e on 
the amount of CDM, the weight of the cosmological con- 
stant, as well as on the weights of the curvatures effects 
encoded on fl r< . and fl a (related to the crossover scale 
and the GB parameter a, respectively). Concerning such 
feature, namely that the phantom like-behaviour in the 
ADGP and ADGP-GB breaks down at the redshifts z b 
and Zb, respectively, can it then be that Zb > 2^? I.e. 
can the GB effect bring the breakdown of the phantom- 
like behaviour in the ADGP to higher redshifts? This 
question has been addressed analytically (in specific sit- 
uations) as well as numerically in more general situations 
and the answer has been always negative (cf. Figs. [6] and 
[7]) : The GB effect seems instead to induce a breakdown of 
the phantom-like behaviour in the ADGP scenario at an 
even smaller redshift. Namely, the ADGP has a regular 
phantom-like behaviour for [0, Zb) whereas the phantom- 
like behaviour in the ADGP-GB model is regular only 
for [0, Zb), with 2% < Zb- Thus, we conclude that the GB 
term does extend the regime of validity of the phantom 
mimicry in the ADGP model. 

We were expecting that new curvature corrections 
would modify the regime of validity of the phantom-like 
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behaviour on the brane. We have shown that this is the 
case by considering a GB term in the bulk. However, 
these bulk curvature corrections rather than enlarging 
the regime of validity of the phantom-like behaviour on 
the brane, they make it smaller. It might be that we 
need to consider exclusively curvature corrections to the 
induced gravity action on the brane, modelled for exam- 
ple through an f(R) term on the brane action [34[. We 
leave this question for a future work. 



Acknowledgements 

The authors are grateful to R. Lazkoz and Y. Shtanov 
for very helpful comments. MBL is supported by the 
Portuguese Agency Fundacao para a Cicncia e Tec- 
nologia through the fellowship SFRH/BPD/26542/2006. 
She also wishes to acknowledge the hospitality of the 
Physics Department of the Universidade da Beira Inte- 
rior during the completion of part of this work. This 
research work was supported by the grant FEDER- 
POCI/P /FIS /57547/2004. 



APPENDIX A: SOLUTIONS OF THE 
FRIEDMANN EQUATION 

The solutions of the Friedmann equation (fl"T|) can be 
written as |35| 



#i = 
H 2 = 

H 3 = 



Si + S 2 - 



1 

3' 



(Al) 



-\(S 1 + S 2 )-^ + i^-(S 1 -S 2 ), (A2) 
-\(Si + S 2 )-^-i^-(S 1 -S 2 ), (A3) 



where 



Si = 



i? + (Q 3 + R 2 



So = 



R-(Q 3 



R l 



(A4) 

Q and R are given in Eq. (|19[) . Only those solutions real 
and positive corresponds to cosmologically interesting so- 
lutions. The latter depends strongly on the ratio between 
the GB parameter a and the crossover scale r c . Once the 
real roots of Eq. (jXTJ) arc identified, it is much more use- 
ful to rewrite them as trigonometric function (if all the 
solutions are real) or as hyperbolic functions if there is a 
unique real solution in order to analyse the cosmological 
evolution of the brane. This is performed in section III. 



APPENDIX B: 



RECOVERY OF THE ADGP 
SCENARIO 



In this appendix, we show that only two of the so- 
lutions of the Friedmann equation of the normal DGP 
branch with GB effects presented in Sect. IHII have a well 



defined limit when a — > 0. These solutions correspond to 
Eqs. ((29j) and {33]). Furthermore, only the solution (|29]) 
reduces to the ADGP scenario when a = 0. 

As it can be noticed, it is not as straightforward to 
take such a limit; i.e. a — ► 0, of the mentioned solutions 
because all the dimensionless parameters involved in the 
problem H, p and b are proportional to a or a 2 . To make 
the task easier, we introduce the following definitions 



p = 




h = 


32 k\ 


(Bl) 


b = 


h a ; 


h = 


8 1 

R' 


(B2) 










(B3) 


H = 


h a ; 


h = 


8 -^H, 
3 r c 



where /1, f 2 and h do not depend on a. 

We start considering the solution (|2"5|) and we make a 
series expansion of tan(#) at a = 0, where 9 is defined in 
Eq. (30), 



tan(6») = 
Therefore, 

6 = 7T 



12/ 1 + 3/|U + 2 (a). 



12/i + 3/ 2 2 )a + 2 (a 



(B4) 



(B5) 



because < 9 < n. 

On the other hand, by combining the last equation and 



cos (| - cr) = i (l + V3c.t) + 2 (x) 
for x — ► and c a constant, we obtain 



(B6) 



cos 



V3 
2 



12A + 3f 2 a 



+ 3 {a). 



Finally, we substitute the last equation in Eq. 
we conclude 



(B7) 
and 



ha 



12A + 3/| a 



+ 




M-l +0 2 (a). 



(B8) 



We equate the lower order in a of the last equation (of 
both handsides of the equation) and we substitute the 
definitions of /1, f 2 an d h- We finally obtain 



H = 




(B9) 



which is the modified Friedmann equation of the non- 
sclf-accclcrating DGP branch that can be obtained from 
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Eqs (1) or (3) (with - sign) with a = and contains the 
ADGP model as a specific solution. 

Following a similar procedure, it can be shown that the 
limit a = of Eq. (03]) 



^=-^( 1 +V 1+ r^j- (B10) 

The limit is well defined although the solution it is not 
physical (at least it is not suitable for the late-time evo- 
lution of the universe). Please notice that Eq. (3) of the 
manuscript with a = is quadratic on the Hubble rates 
and therefore has two roots. This second root; i.e. Eq. 
(BIO), after a time reversal {t by —t) corresponds to the 
self-accelerating DGP solution. 

On the other hand, the solution (32) of the paper has 
not a well defined a = limit because the left hand side of 
the equation is proportional to a while the right handsidc 
lower order in a is a constant different from zero. 

If one tries to take the limit a — > of the solution 
presented in Eq. (25) it turns out that r\ must be complex. 
In fact, in that case one is getting an analytical extension 
of the solution which corresponds to the solution ([29]) . 
The reason behind this behaviour is that the solutions 



(j29|) and (25) are related by a sort of "Wick rotation" 
where 6 = ir/. Notice at this respect that the function N 
changes sign and the prolongation is well defined. 

For the rest of solutions presented in section III, al- 
though all of them comes from the general solution (Al), 
the limit a — > cannot be taken because 1/4 < b. 

In summary, the recovery of the ADGP scenario from 
the normal DGP branch with GB effects requires that 
(i) the value of the parameter b, which is proportional 
to a, to be such that < b < 1/4 and (ii) the dimcn- 
sionlcss amount of the branc total energy density has to 
fullfill p < pi . Therefore, in order to recover the ADGP 
scenario from the model we propose, the GB parameter 
has to be bounded from above and the brane total energy 
density has also to be bounded by a maximum threshold. 
Along section IV, we have restricted to this solution as it 
is the most favourable solution obscrvationally because 
the parameter b is expect to be small which implies that 
the most likely set of cosmological parameters that char- 
acterises the model would lie on the bottom of the left 
handside corner of figure 2. This last feature on itself im- 
plies that the brane energy density cannot be too large, 
more precisely p m i n < p\ which implies in particular that 
there is a range of values of p such that p\. 
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